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The conventional thermal instability criterion can not be applied to the advection- 
dominated accretion disks around black holes where the radiative cooling is insuffi- 
cient to balance the viscous heating. The surface density change associated with the 
temperature perturbations, which was usually neglected in deriving the conventional 
criterion, was recently shown to be much significant in the advection- dominated 
disks. Considering both advection and surface density change, I suggested a gener- 
alized thermal instability criterion. By applying it to the optically thin and optically 
thick advection- dominated disks, I found that the former one is thermally stable 
and the latter one is thermally unstable against short wavelength perturbations, 
which agrees well to those found recently by both analytic and quantitative stability 
^ ! analyses. 
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Q\ . Accretion disk around black holes can provide a large amounts of released energy and is 

thus believed to be the source of high energy emission of Galactic black hole X-ray binaries and 
active galactic nuclei . The stability properties of accretion disks are very important because 
a global violently unstable disk may not exist in nature and some instabilities restricted in 
a certain region of the disk may contribute to the observed light variations in many active 
astronomical systems. Shortly after the construction of the standard model of an optically 
thick, geometrically thin accretion disk 2 , the inner region of such a disk was found to be both 
thermally and secularly unstable 3,4 . The hot optically thin disk model, proposed by Shapiro, 
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Lightman & Eardley to explain the hard X-ray spectra of black hole X-ray binary Cyg X-l, 
was also proved to be thermally unstable 6 . The general stability criteria of accretion disks 
were derived by Piran 7 who considered various kinds of viscosity laws and cooling processes. 
These criteria can be used as an a priori check if some choice of cooling mechanisms gives rise 
to a stable disk even before the complete disk structure equations are solved. The thermal 
instability criterion can be expressed as 8 : 

dlnQ- dlnQ + 

{ -dMT h < ( ^T" )E ' (1) 

where Q + , Q~ and T are the viscous heating rate, radiative cooling rate and temperature. 
£ is the surface density defined as £ = pH where p and H are the density and disk height 
respectively. Such a criterion is identical to that appeared in Piran 7 if we relace T with 
H . The thermal instability can be well understood according to this criterion. Considering a 
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small temperature increase and assuming that the surface density does not change within the 
thermal time scale, if the viscous heating rate increases more rapidly than the radiative cooling 
rate does (as indicated in inequality (1)), it will lead to a further increase of temperature and 
the thermal instability sets in. If the standard viscosity law is adopted, Piran has pointed out 
that the inner region of accretion disks is always thermally unstable no matter what kind of 
cooling mechanism is involved 7 . 

However, the thermal instability criterion we mentioned above can not be applied directly 
to the advection-dominated disks. In the advection-dominated disk models, the radiative 
cooling is inefficient and most of the energy generated by viscous heating is transfered to 
increase the entropy of accretion flow and then is advected inwards onto the black hole 9 > 10 > u . 
To derive a criterion for the thermal instability of advection-dominated disks, we should 
consider the advected energy in the energy conversation equation. Moreover, the surface 
density change associated with the temperature perturbations, which used to be negligible in 
previous stability analyses of the geometrically thin disks, should be considered in deriving a 
thermal instability criterion of advection-dominated disks. Some recent works on the stability 
of advection-dominated disks have indicated that the surface density change is much significant 
in these disks comparing to the minor pressure change 12 > 13 > 14 . In this letter, I will present 
a generalized thermal instability criterion by considering both advection ad surface density 
change. 

In deriving the conventional thermal instability criterion, the energy balance between heat- 
ing and cooling was assumed in the unperturbed case. However, in an advection-dominated 
disk, this balance is realized among viscous heating, radiative cooling and advection (for sim- 
plicity, the effects of magnetic field and thermal diffusion are not included in the present 
study). That is, 

Q+ = Q' + Q ad (2) 

where Q + , Q~ and Q ad are the viscous heating rate, radiative cooling rate and the energy 
advection rate respectively. The crucial point of the thermal stability of accretion disks is to 
examine whether or not this balance is still held if the temperature perturbations are involved. 
To do this examination, we first write some relative energy changes as follows: 

5Q + 5T 5Q- 5T 5Q ad 5T 

where the relative changes of Q + , Q~ and Q ad are expressed as the functions of temperature 
and surface density. Xi and yi (i = 1, 2, 3) can be defined as: 

dlnQ + dlnQ + dlnQ~ dlnQ~ 

MnQ ad ^ ,dlnQ ad , 
X3 = ( ^^ )s ' m = i ^Z )T (4) 
If we also define q = Q ad /Q + and Q~ = (1 — q)Q + , we can obtain: 

SQ+ _ S Q- _ 5Q ad 

— = {( Xl _ X3 ) - ( 1 - q ) ( X2 - X3 )\ — + [( yi - y 3 ) - (1 - q) (y 2 - y 3 )] — (5) 
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The thermal instability will arise if the left side of above equation is positive, which means 
that with a small increase of temperature the viscous heating rate grows more fast than the 
sum of radiative cooling rate and energy advection rate does. Therefore, a thermal instability 
criterion including the advection term and the surface density change will be: 

xx - (1 - q)x 2 - qx 3 > -8tt[v\ ~ (1 - q)y2 ~ qy?] (6) 

where 5st = ^/^i which describes the relative changes of surface density to the tempera- 
ture perturbations. This criterion can be regarded as the generalization of the conventional 
criterion (inequality (1)). If we keep the surface density as a constant and take q as zero, it 
will be the identical to the inequality (1). Because the advection term and the surface density 
change are considered in our new criterion, we can use it now to check the thermal instability 
of advection-dominated disks. 

For an optically thin disk, the pressure is dominated by gas and the opacity is dominated by 
electron scattering 10 . We assume that radiative cooling is provided by thermal bremsstrahlung 
with emissivity as 

Q~ = 1.24 x 10 21 Hp^^ergss^cm- 2 . (7) 
The viscous heating rate is given by the standard formula 

Q+ = S,(r§) 2 , (8) 

Where v is the viscosity parameter and Q is the angular velocity. The advection cooling rate 
is taken as 

Q ad = C v [Y,V r °— - (r 3 - ^)TV r °—l (9) 

where C v = 3c 2 /2T, T 3 = 5/3, V r is the radial velocity and c s is the local sound speed. 
Considering p = KpT, H = c s /Qk, c 2 — T>l Pi v = olc s H and V r ~ u/r, we have Q + oc ET, 
Q- oc S 2 and Q ad oc ST 2 [s(T) - s(E)]. Here s(T) = ^ and s(E) = and we assume 
them to be constant, which is valid if T and E are the power law functions of r as indicated, 
for example, by the self similar solutions 9 . Therefore, for an optically thin disk we have 
x\ — 1, yi — 1, X2 = 0, i/2 = 2, X3 = 2 and y 3 = 1. Taking these values into inequality (6), we 
can obtain a criterion for thermal instability of an optically thin disk, that is 

5st < 1 - zr—. (10) 

1 - q 

If the optically thin disk is radiative cooling dominated, which means q ~ 0, above in- 
equality requires 5^t < 1- In the long perturbation wavelength case, there is no appreciable 
surface density change so we have 8y,t ~ 0. In the short perturbation wavelength case, a recent 
study by Wu 14 showed —1 < 5^t < 0. Therefore, our criterion proves again that an optically 
thin disk is thermally unstable if it is radiative cooling dominated. On the other hand, if 
the optically thin disk is advection-dominated, which means q ~ 1, inequality (10) requires 
<5et < —00. Since we always have — 1 < 5y,t < for an optically thin, advection-dominated 
disk, obviously this inequality means that the disk is thermally stable. 

If the optically thin disk is two-temperature one, the radiative cooling may be dominated 
by Comptonization through the loss of energy of electrons. Then the radiative cooling takes 
the form of 5 : 

Q_ = (4kT e /m e c 2 )pHK es U r c. (11) 
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where m e , T e are the mass and temperature of electron, K es = 0Acm 2 g~ 1 is the electron 
scattering opacity. U r is the radiation energy density of soft photons, which we assume, for 
simplicity, does not change on time scale shorter compared to fi -1 . Because the ions and 
electrons are coupled by collisional energy exchange, the loss of energy of electrons can be 
balanced by the energy capture from ions. The exchange rate is 

Qi- = {'S/2) P Hv E k{T i - T e )/m p (12) 

where v E is the electron-ion coupling rate, and can be approximated by v E = 2 Ax 10 21 ZnApT e ~ 3 / 2 
where the Coulomb logarithm InA is about 15. The ion temperature T« is usually one or two 
orders higher than the electron temperature T e . Taking Q_ = Qi_ and T« >> T e , we can get 
Q_ oc E 7 / 5 ^ 1 / 5 . Then we have x 2 — 1/5 and y 2 = 7/5, then the thermal instability criterion 
now becomes: 

^ <2 -W^y < 13 > 

Similar as we discussed above, this inequality clearly shows that a radiative cooling dominated 
two-temperature disk is thermally unstable and an advection-dominated two-temperature disk 
is thermally stable. A more recent quantitative stability analysis also reaches the same con- 
clusion 15 . Moreover, from inequalities (10) and (13), we note that in the long wavelength 
perturbation case where 5y,t ~ 0, the bremsstrahlung cooling disk and the two temperature 
disk may become thermally stable when q > | and q > |, respectively. 

For an optically thick accretion disk, if it is dominated by radiation pressure and electron 
scattering, we can write the viscous heating rate as Q + = Sz/(r^) 2 cx and the radiative 
cooling rate as Q~ = 4 °^ 4 oc where k is the opacity dominated by electron scattering. 
The advection cooling rate, which is described by equation (9), can be expressed as Q ad cx 
^3-[s(T) — s(E)] for an optically thick disk. Therefore, we have x-y — 8, y± — —1, x 2 — 4, 
y 2 = — 1, x% = 16 and y% = —3. The criterion for thermal instability of optically thick disks 
can be written as: 

2 

5 ST > 6 - -. (14) 
Q 

If the otically thick disk is radiative cooling dominated, above inequality becomes 5^t > 
—00, which is satisfied without any doubt. Therefore, this indicates that the optically thick, 
radiation pressure and radiative cooling dominated disk is thermally unstable, which was first 
proved by a linear study of Shakura & Sunyaev 3 . If the optically thick disk is advection- 
dominated, the thermal instability criterion requires 5^t > 4. A recent study of Wu 14 
pointed out that for an optically thick, advection-dominated disk we have <5et ~ 8 in the 
short wavelength perturbation case, which is the result of the fact that there is no appreciable 
vertical integrated pressure change associated with the temperature perturbations. Therefore 
we still obtain that the optically thick, advection-dominated disk is thermally unstable, even 
if the effects of thermal diffusion and magnetic fields are not included. 

In addition, we note that in the long wavelength perturbation limit where we have 5^t ~ 0, 
the optically thick disk may become thermally stable when q > |. This also agrees well with 
the implication of the turning point from the radiative cooling dominated branch to the 
advection-dominated branch in the thermal equilibrium curve 10 . 

The generalized criterion can be applied to investigate the thermal stability properties of 
both radiative cooling dominated disks and advection-dominated disks. The results qualita- 
tively agree well to those obtained by some recent detailed stability analyses, which shows 
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that an optically thin advection-dominated disk is thermally stable while an optically thick 
advection-dominated disk is thermally unstable against the short wavelength perturbations 
12 ' 13 . Moreover, without solving the detailed disk structure, our new criterion can be used to 
check whether ot not the disk is thermally stable when the advection effect and the different 
cooling mechanism are involved. 

We note that our new criterion depends on the relative surface density change <5et, which 
is a function of the perturbation wavelength and other disk parameters such as viscosity, Mach 
number, disk height and advection term 14 . In the limit of short or long wavelength perturba- 
tions, we can easily obtain the thermal stability properties of the disks from the criterion by 
adopting the known values of 5st- However, in the case of intermediate perturbation wave- 
length, we have to obtain first by numerically solving the perturbed equations before we 
adopt the new criterion. In this sense, we think that the generalized thermal instability crite- 
rion suggested in this letter is only appropriate in the short and long perturbation wavelength 
limits and the quantitative stability analyses of accretion disks are still needed for the case of 
intermediate perturbation wavelength. 

I thank Dr. Xingming Chen at Lick Observatory for helpful discussions. 
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